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Abstract 

Q. 

We investigate a simplified version of the ACOT prescription for calculating deeply 
inelastic scattering from Q 2 values near the squared mass Mfj of a heavy quark to Q 2 
■ much larger than Mjj. 
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1 Introduction 



The inclusion of heavy quark effects in deeply inelastic scattering (DIS) is an interesting 
theoretical problem involving two hard scales in a perturbative analysis. This issue is also 
of phenomenological importance. The charm contribution to the total structure function F 2 
at small x, at HERA, is sizeable, up to 25%. Through the charm contribution to scaling 
violations, the treatment of charm also has a significant impact on the interpretation of fixed 
target deeply inelastic scattering data. Thus, a proper description of charm contributions to 
deeply inelastic scattering is required for a global analysis of structure function data and a 
precise extraction of the parton densities in the proton. 

At scales Q <> Mh, the contribution to deeply inelastic scattering of a heavy quark of 
mass Mh can be calculated in the so-called fixed-flavor-number (FFN) prescription from 
hard processes initiated by light quarks (u, d,s, . . .) and gluons, where all effects of the 
heavy quark (H) are contained in the perturbative coefficient functions. This prescription 
incorporates the correct threshold behavior, but for large scales, Q ^> Mh, the coefficient 
functions at higher orders in a s contain potentially large logarithms In 1 (Q 2 / M H ) , which may 
need to be summed 0, ^, Such a summation can be achieved by including the heavy 
quark as an active parton in the proton. The simplest approach incorporating this idea is 
the so-called zero-mass variable-flavor-number (ZM-VFN) prescription, where heavy quarks 
are omitted entirely below some scale Q ~ M H and included as massless partons above 
this threshold. This prescription has been used in global analyses of parton distributions for 
many years, but it has an error of order M H /Q 2 and is not suited for quantitative analyses 
unless Q > M H . 

Considerable effort has been devoted to including heavy quark effects in deeply inelastic 
scattering in such a way that the calculated structure functions match those of the FFN 
prescription in the region Q « Mh while they match those of the ZM-VFN prescription 
for Q ^> Mh- Two prescriptions of this sort, the Aivazis-Collins-Olness-Tung (ACOT) 
@] and the Thorne-Roberts || prescriptions have been used in recent global analyses of 
parton distributions [EL f7|. More recently, additional variable-flavor- number prescriptions 
with non-zero mass have been defined in the literature |§. If one could sum perturbation 
theory, the calculated structure functions should be identical for any prescription that does 
not neglect the heavy quark mass. However, the way of ordering the perturbative expansion 
is not unique, so that the results generally differ at any finite order in perturbation theory. 

In this paper, we will investigate a modification of the ACOT prescription advocated by 
Collins ||. It has the advantage of being simple to state and of allowing relatively simple 
calculations. This simplicity should be convenient for phenomenological analyses at the Born 
level. In addition, it could be crucial for implementing a variable-flavor-number prescription 
with non-zero mass at next-to-leading order in global analyses of parton distributions. 
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2 Schemes and partons 



We consider the structure function F 2 (x, Q). The observable structure function can be writ- 
ten in terms of parton distribution functions / and a calculable partonic structure function 
F 2 as 

F 2 (x, Q) = Cdi £ / a/p (£, /i) F 2 (a, x/(, fi/Q, a s (fi)) + 0(A 2 /Q 2 ). (1) 

Jx a 

This factorization formula has corrections of order A 2 /Q 2 , where A is a typical scale of 
hadronic physics, perhaps the mass of the rho meson. Here \x is the renormalizat ion/fac- 
torization scale. (For simplicity, we do not distinguish these two scales.) The partonic 
F 2 has a perturbative expansion in powers of a s (fi). There is a sum over parton types 
a — g,u, u, d,d, . . ., and there is a function f a / p for each parton. 

To begin, we must pick a scheme for the definition of a s and for the definition of the 
parton distribution functions. In fact, we will use multiple schemes, following the prescription 
of ||. Each scheme is designated by the number N of "active" quark flavors. Scheme N is 
designed to be most useful for physical scales Q in the range Mjy ^ Q ^ Mjv+i, where M/v is 
the mass of the iVth quark flavor. 

We define what we mean by the running coupling (//) in scheme N by defining how 
we perform renormalization. We renormalize using the CWZ prescription |TIJ. Briefly, 
divergences involving active parton loops are removed with an MS subtraction, but when 
active parton external lines couple to a loop containing "non-partonic" lines (quarks iV + 
1, N+2, . . .) the renormalization is by subtraction at zero external momentum. The running 
of (//) is controlled by the usual MS renormalization group equation with the contributions 
from quarks 1, . . . , N in the beta function. 

In the scheme iV there are parton distributions for gluons and for quark flavors 1, . . . , N, 
but not for quark flavors N+l, .... The parton distribution functions are defined to be proton 
matrix elements of certain operators ||11|| . The operator products are ultraviolet divergent 
and are renormalized according to the CWZ prescription. They obey a renormalization 
group equation - the usual DGLAP equation - in which contributions from quark flavors 
N + 1, . . . do not appear in the kernel. 

Note that the parton distribution functions are non-perturbative objects. There is no 
question of neglecting any masses in the definition. The quark masses do not appear in the 
kernel of the evolution equation, but this is because renormalization counter terms are mass 
independent, not because the parton distributions themselves are mass independent. 

Note also that by using CWZ renormalization throughout, quark loops for the non- 
partonic flavors decouple from calculations when the momentum scale is small compared to 
M N+ i. Thus for Q <C M N+ i it is a good approximation to leave the non-partonic flavors out 
of calculations altogether. (But leaving heavy flavors out of calculations is not part of the 
definition of scheme N as used here; it is a separate approximation.) When one does leave 
the non-partonic flavors out of calculations, the remaining renormalizations are via the MS 
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prescription. Thus one commonly refers to parton distributions in the prescription described 
here as MS parton distributions. 

There is a perturbative connection between schemes N and N+l which we can represent 

by: 

CW) = f? /p {x,n) + £ f l ^A ab (x/{;,»/M N+1 ,a s (ri) (2) 

b Jx $ 

Here the index a runs over g, the N light quarks and their anti-quarks, and the heavy quark 
H and anti-quark H. For a — H, H we define fa/ p { x i^) = in the first term on the right 
hand side. The index b runs over only the gluon and light quarks. At order a s , the only 
nonvanishing Am, is Affg- Thus we may say that the H distribution arises perturbatively 
from g — ► HH. 

The a* terms in the perturbative expansion of the A ab vanish at fi — M^+i- Thus one 
derives the simple matching conditon jn|: Mjv+i) = fa/pi x i Mat+i) + C(«s). The 

analogous connection is known at order but we do not repeat it here.Q 

3 Why is there an ambiguity? 

When we construct the hard scattering cross section in the presence of heavy quarks, there 
is a factorization ambiguity that is not present in the case of light quarks. Let us see why. 
Consider two enormously simplified examples that illustrate the principle. 

First, suppose that there are only gluons and one light quark L with mass M^A. 
Suppose, in addition, that the light quark is its own anti-particle, so that a quark L and the 
anti-quark L are the same particle. Both the gluon and the light quark are considered to be 
active partons. In order to simplify the notation, let us take a moment J dxx n F 2 (x,Q) of 
F 2 , so that we get a factorization formula involving the corresponding moments of F 2 and 
of the parton distributions. The dependence on the moment number will not be indicated. 
To further simplify the notation, let us set the factorization and renormalization scale fi to 
Q. Then 

F 2 (Q) ~ F 2 {L, a s (Q)) f L/p (Q) + F 2 {g, a s (Q)) f g/p {Q). (3) 

Now F 2 is an observable, so its definition is fixed. We have defined the parton distribution 
functions, and the two parton distribution functions are independent. Thus this equation 
all but fixes the definition of F 2 (L, a s (Q)) and F 2 (g, a s (Q)). The only possible modification 
would be to add terms proportional to powers of Ml/Q 2 to F 2 (a, a s (Q)), at the cost of 
subtracting the same terms from the power suppressed remainder, 0(A 2 /Q 2 ), in the factor- 
ization formula (Eq. ([[])). The simplest solution, which is uniformly adopted, is not to allow 
a Ml/Q 2 dependence in F 2 (a, a s (Q)). 

Now suppose that there are only gluons and one (self-conjugate) heavy quark H with 
mass Mn ^> A. Then in the scheme in which both the gluon and the heavy quark are 
1 At order a 2 s , the matching at /i = M^+i is no longer continous in general 
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considered to be active partons we have 

F 2 (Q) ~ F 2 (H, M H /Q, a s (Q)) f H/p (Q) + F 2 (g, M H /Q, a s (Q)) f g/p (Q). (4) 

Here the F 2 (a, Mh/Q, a s (Q)) (a=g,H) depend on the heavy quark mass and we cannot move 
terms of order M H /Q 2 into the power suppressed corrections because only terms of order 
A 2 /Q 2 are allowed there. Thus it seems that we have no freedom. But we do: fu/p{Q) and 
f g /p(Q) are not independent. Since, according to Eq. (g), heavy quarks evolve from gluons, 
we have a relation of the form 

f H / P (Q) = V H/g (ln(Q/M H ), a s (Q)) f g/p {Q). (5) 

Here V has a perturbative expansion that is obtained by solving the evolution equation. The 
first term has the form V ~ a s 7 ln(Q/M) where 7 is a constant. Using Eq. fl5|) in Eq. @ 
we obtain 

F 2 (Q) ~ {P 2 (H, M h /Q, a s (Q)) V H/g (HQ/M H ), a s (Q)) + F 2 (g, M H /Q, a s (Q))} f g/p (Q). 

(6) 

Evidently, there is some freedom to move pieces from the first term in braces to the second. 
There is a constraint. For Mh/Q — » 0, it is possible to neglect Mh in the calculation of 
the F 2 . On the other hand, V(\ii(Q/Mh), a s (Q)) does not have a smooth Mh/Q — > limit. 
This is not a problem in applications because by solving the evolution equation one sums 
the leading logarithms [a s ln(Q/M#)] n in V. It is important not to undo this summation. 
Thus the factorization scheme that we adopt should have the property that the functions 
F 2 (a, Mh / Q,a s {Q)) have a finite limit as Mh/Q — > 0. This still leaves us the option of 
adding a term like c x (M H /Q 2 ) to F 2 (H,M H /Q,a s (Q)) and subtracting c x (M H /Q 2 ) V 
from F 2 (g,M H /Q,a s (Q)). 

Suppose now that we have calculated F 2 (g,M H /Q,<y s (Q)) and F 2 (H,M H /Q,a s (Q)) in 
some convenient prescription - for example the prescription analyzed in || based on "on- 
shell" heavy quarks. Then we could define a new prescription with 

Fr(H, Mh/Q, a s (Q)) = F° ld (H, 0, a s (Q)) (7) 

and 

Pr(9, Mh/Q, a s (Q)) = F° ld (g, M H /Q, a 8 (Q)) 

+ {F° ld (H,M H /Q,a s (Q)) - F° ld (H,0,a s (Q))} x V H/g (ln(Q/M H ),a s (Q)). (8) 
In the following subsection, we shall give a prescription based on this observation. 

4 A prescription for resolving the ambiguity 

Having seen the main idea, let us put the parton indices and the momentum fraction variables 
back. Suppose that there are N + 1 quark flavors that we consider to be active. Let the 
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heaviest active quark, quark N + 1, be labelled H. Suppose that quark H has a mass that is 
large compared to the hadronic mass scale A. In the N + I flavor scheme the factorization 
equation is 

F 2 (x,Q) ~ / £/«/p(£^) F 2 {a,x/£,M H /Q,n/Q,CL a <ji)), (9) 

where the sum over a includes a = H and a = H. 

The parton distributions for the gluon and iV + 1 quark flavors are obtained from the 
distributions in the scheme with only N quark flavors. First we use the perturbative matching 
relation (Eq. (||)) at a scale near /i = M H , then we use the evolution equations to give the 
distribution functions at scale \i. Since there are more output functions than input, we 
obtain a perturbative relation giving the heavy quark distribution functions at scale \i in 
terms of the light quark and gluon distribution functions at the same scale. This relation 
has the form 

rl df 

fH/p{t,n)= — £ fa/ P (r,fi) V H/a {£/r;fi/M H ,a s {fi)), (10) 

J £ T I T T T~T 



a^HH 



with an analogous equation for the heavy anti-quark H. Inserting this relation into Eq. (|S|). 
we have 

F 2 {x,Q)~f\ J2 fa/ P {L»)T a {x/Z,M H /Q,fi/Q,a s {fi)), (11) 



a^H,H 



where 



T a (z,M H /Q,/j,/Q,a s (iJ,)) = F 2 (a, z, M H / Q, n/Q, a s (/j,)) 

+ j dX F 2 {H, z/X, M H /Q, fx/Q, a s {fi)) V H/a {X, fx/M H , a s (fx)) 

+ J dX F 2 {H, z/X, M H /Q, /i/Q, a a (n)) V H/a {X, fx/M H , a s {fi)). (12) 

We see that we have the same situation as in the simple example given earlier. Taking \x 
to be of order Q, one can shift contributions of order M H /Q 2 between the hard scattering 
functions F 2 for H and H and the corresponding functions for the light quarks and the gluon 
while keeping the functions T unchanged for each light quark and gluon flavor a and without 
ruining the property that all of the functions F 2 have finite Mh /Q — » limits. 

This freedom can be exploited to make the calculation of the functions F 2 simpler. In 
particular, we can adopt a prescription proposed by Collins 0: 

Simplified ACOT (S-ACOT) prescription. Set M H to zero in the calculation of 
the hard scattering functions F 2 for incoming heavy quarks. 

This observation tremendously simplifies the calculation of F 2 for a = H as it reduces to 
that of the light-quark result. For example, at order a]., the F 2 functions for heavy quarks 
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and light quarks are independent of Mh, and the calculation reduces to that of the simple 
massless result. The a] gluon contribution to F 2 acquires an Mh dependence when the gluon 
couples to a heavy quark loop, which is probed by the virtual photon (g / y —> HH). 

Note that the hard scattering functions F 2 obey a renormalization group equation that 
is different from the standard renormalization group equation that applies in the case of 
light flavors only. To see why this is so, imagine that the F 2 functions obeyed the usual 
renormalization equation and that F 2 (H, zj A, M H /Q, /i/Q, a s (A0) were independent of M H 
at some fixed value of fi. Then F 2 for a = H,H would depend on Mh at other values of \x 
because the standard renormalization group equation mixes the heavy and light flavors and 
would mix M H dependence into F 2 (H, zj A, M H /Q, [i/Q,a s (/i)). 

This Simplified ACOT prescription has the advantage of being simple to state. In ad- 
dition, its calculational simplicity could be crucial for applying variable-flavor-number pre- 
scriptions with mass in analyses that go beyond first order in a s . 



5 The ACOT and S-ACOT prescriptions at first order 

In this section, we analyze the ACOT prescription and its simplified version, the S-ACOT 
prescription, at order a]. We consider neutral current deeply inelastic scattering from a 
proton target in a Q 2 regime in which quarks 1, . . . , N can be considered as light while a 
single quark, H, is considered to be heavy or light depending on the value of Q 2 . In order 
to keep the analysis simple, contributions from all heavier quarks are ignored. We further 
simplify the problem by supposing that the vector boson current that probes the proton 
couples only to the heavy quark H and its anti-quark H, but not to the lighter quarks or 
gluons. (In a crude approximation, this is like taking ordinary deeply inelastic scattering but 
demanding that H appear in the final state.) Let F be one of the structure functions F\, 
F 2 /x or F 3 for our special vector boson. We choose the factorization and renormalization 
scales n equal to Q. We write the factorization formula for F in the shorthand notation 

F = J2F a ®f a/p + 0(A 2 /Q 2 ), (13) 

a 

where the sum runs over all partons a including H and H. The ® denotes a convolution, so 
that Eq. ( ji3|) means 

F(x, Q) = E / f ^Wt> Q > a ^ fa '^ + °^ 2 /Q 2 )- (14) 

The functions F a in Eq. (|13|) are the hard scattering functions. They have an expansion in 
powers of a s . Keeping the first two terms in this expansion gives 

F = F^®f H / P + Ff®f H/p 

+ ® f B / P + F ( h ] ® f H/p + F H l) ® f H/p 

+ 0(a 2 )+0(A 2 /Q 2 ), (15) 
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where F^ is the order a™ contribution to F a . 

In the ACOT prescription, we choose the order zero hard scattering function for a heavy 
quark to be 

FS\coT = FS\M H ), (16) 

where F^ (Mh) is the calculated structure function for scattering from an initial state heavy 
quark of mass Mh that is on its mass-shell. (To be precise, the heavy quark transverse 
momentum is taken to be zero, and the structure functions Fi, F 2 /x, F 3 are extracted from 
the tensor W^ u using the usual formula but with the quark momentum k^ replaced by a 
light-like vector k^ = k^ — [M H /(2u ■ k)] u^, where is a light-like reference vector in the 
plane of g M and k^ 1 .) The function F^(Mh) is rather complicated (see Ref. ||, [13], [14|])> so 
we do not reproduce it here. The definition of Fft\coT for a heavy anti-quark is analogous. 
The order a s gluon hard scattering function in Eq. (|15|) has three pieces: 

^acot = FU(M H ) - F$\M H ) ® fSj g {M H ) - Fj°\M H ) ® fff /g (M H ). (17) 

Here F^ (Mh) is the calculated structure function for scattering from an initial state massless 
gluon that is on its mass-shell, using graphs with a heavy quark loop and a heavy anti-quark 
loop. The function f^J g (Mn) is the calculated a] contribution to the distribution of heavy 
quarks in an on- mass-shell gluon. Similarly, f^ g (M H ) is the calculated a\ contribution to 
the distribution of heavy anti-quarks in an on-mass-shell gluon. Both functions are given by 

f§) B &Q) = f§tit,Q) = ln |r P ^ (18) 

where P q / g is the usual gluon — > quark splitting function P qg (£) = Tp (£ 2 + (1 — £) 2 ). 

The first order hard scattering function for a heavy quark has a structure similar to that 
of the corresponding function for a gluon, 

^acot = F { h\Mh) - FP{M H ) ® fSj H (M H ). (19) 

Here F h ~\Mh) is the order a] contribution to the structure function for scattering from an 
initial state heavy quark that is on its mass-shell, as given in fl3| . The function f^J H (Mn) is 
the calculated a), contribution to the distribution of heavy quarks in an on-mass-shell heavy 
quark, 



ir r~\\ n a s(Q) 



(20) 



2tt 

where the + subscript denotes the usual prescription, 

[>(£)] = "/(I)} fXO- (21) 

Jo 1 J + Jo 

This result can be calculated easily from the definition [TT| of MS parton distribution func- 
tions or it can be extracted from the ACOT subtraction terms in | 13fl . 

This defines the ACOT prescription at order a]. The prescription has two important 
properties. 
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p(o) 

r H,ZM-VFN 


= 4 0) (o), 


r g,ZM-VFN 




r H,ZM-VFN 


= 4 1} (o) 



• Property 1. For Q/M — > oo, the hard scattering functions in Eq. (|T^) approach the 
hard scattering functions of the ZM-VFN prescription, in which the heavy quark H is 
taken parton with zero mass. 

To be specific, we have the relations: 



4 0) (o) ® fSj B (o) - 4 0) (o) ® /^(o), 

4>) (22) 

We should note that in a calculation with M# = 0, there are infrared divergences. Thus the 
calculations are performed in 4 — 2e dimensions of space-time and 1/e poles appear. The 
pole terms cancel in Eq. fl22"[) . Instead of using dimensional regulation and taking e — > 0, 
one could use an infrared regulator mass m and let m — > 0. 

• Property 2. For Q of order Mg, the structure function is that of the fixed- flavor- number 
prescription, up to corrections of order a 2 . 

In the fixed-flavor-number prescription we have 

F = FW(M H ) ® / a/j> + 0(a 2 s ) + 0(A 2 /Q 2 ). (23) 

Property 2 follows because when the factorization scale, \i — Q, is of order M, the heavy 
quark distribution function is given by the perturbative formula 

f H , P = f ( H} a ®f 9 /p + 0{a 2 s ). (24) 

First of all, since fjj/ p is of order a s , the terms 

F { H ] ®f H , P + F { ff®f R/p (25) 

in Eq. ( |T5| ) are of order a 2 and can be dropped. This leaves 

F = Fg\M H ) ® / H/P + 4° } ( M ^) ® /h/p 

+ [^ (1) (Mh) - Ff{M B ) ® fS) g (M B ) - Ff{M H ) ® /$,(M ff )] ® / a/p 

+ 0(a 2 ) + 0(A 2 /Q 2 ). (26) 



Inserting Eq. (p4|) into Eq. (pq), we obtain Eq. (pq). 

We can summarize properties 1 and 2 by saying that the ACOT prescription interpolates 
between the zero-mass prescription for M B -C Q and the fixed-flavor-number prescription 
for M H ~ Q. 
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What about the simplified ACOT prescription? We use the same formulas as for the 
ACOT prescription, but set Mh = in the hard scattering functions for heavy quarks: 

p(o) _ p(0)/ n \ 

-^S-ACOT — r H l u J> 

^S-acot = F^\M H ) - 4 0) (0) ® fSj g (M H ) - 4 0) (0) ® fjp /g {M H ), 
^s-acot = 4 1} (0) - 4 0) (0) ® fS} H (0). (27) 

Repeating the derivation just given, we see that properties 1 and 2 hold for the S-ACOT 
prescription. That is, the S-ACOT prescription also interpolates between the zero-mass 
prescription for Mh <C Q and the fixed-flavor-number prescription for Mh ~ Q- However 
the S-ACOT prescription has the advantage that the functions Fjj are easier to calculate. 



6 Comparison of different prescriptions 

In this section, we investigate how well the matching properties among the different prescrip- 
tions work. In the plots presented, the masses and couplings have been fixed to be consistent 
with the values used in the CTEQ4L/M fits, i.e. m c = 1.6 GeV, m;, = 5 GeV Below 
the bottom threshold nf = 4 active flavors are used for a s , and the scale has been chosen as 

fi = Q- 



6.1 Parton distribution matching 

In the threshold region, the structure function calculated in the S-ACOT prescription matches 
that calculated in the FFN prescription, in which the heavy quark does not appear as a par- 
ton. This matching, Property 2, results from the fact that the heavy quark distribution 
function fH/ P (x,fi), matches the approximate function 

Ih/M = f ( H/ 9 (»)®f 9 /p(»)- (28) 

Both functions vanish at \i — Mh and, for fi > Mh, the difference between them is of order 
a 2 s . In this subsection, we study the quality of this matching. 

In Figs. 1-3 we plot fs/p( x , /-0 an d fii/ P (x, A*) for the case of the charm quark. In Fig. 1, we 
use CTEQ4L parton distributions, which are based on lowest order evolution. Fig. 1 reveals 
that the matching works very well when the order a] evolution kernel is used for /h/ p (x, //) 
and the order a l s perturbative expression (Eq. (|2"%D) is used for /. In Fig. 2, we use CTEQ4M 
parton distributions, which are based on NLO evolution. We observe a mismatch because 
the full NLO evolution kernel is used for fn/ p while the order a] perturbative expansion is 
used for /. The difference between / and / is of order a 2 s) but this difference is numerically 
quite large. Finally, we see in Fig. 3 that a close match is restored when the NLO evolution 
kernel is used for / while / is defined using a modified version of Eq. (p?f ) in which we replace 
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Figure 1: CTEQ4L charm quark density fii/p{x,Q) and the approximate form fH/p(%,Q), 
Eq. fl28|), at x = 0.05 as function of Q. Since Eq. (^8]) is an order a\ perturbative approxi- 
mation and the evolution kernel for the CTEQ4L parton distribution is also order a\, these 
curves match closely at threshold. 

the leading order g — > H evolution kernel in Eq. (|i~8| ) by the next-to-leading order g — > H 
evolution kernel.^ 

We can draw two conclusions. First, the threshold matching discussed in the previous 
section will work order by order for the perturbation expansion.^ Second, the order a 2 s terms 
in the evolution kernel are quite large, so that the leading order calculations illustrated in 
this paper may not be sufficient for obtaining accurate predictions.^ 

6.2 Structure function matching 

In this subsection, we examine predictions for F 2 c (x,(5), which we define here to be the 
contribution to F 2 (x,Q) from graphs in which the current couples to a charm quark. We 
compare F£(x, Q) calculated with the S-ACOT prescription at order a s with that calculated 
with the original ACOT prescription, the ZM-VFN prescription in which the charm quark 
can appear as a parton but has zero mass, and the FFN prescription in which the charm 

2 For the NLO perturbative expansion, we have included the a 2 splitting kernels, For simplicity, 

we have ignored iterated terms such as P^jl which contribute as hi 2 (/i/Mh ), and hence only play a 

role away from threshold. 

3 At order a 2 , Property 2 is still preserved even though the matching conditions on the parton distribution 
functions are modified. The a 2 s matching conditions shift the evolved parton distribution functions fa /pi 
but they also shift the "perturbative" parton distribution functions fn/p leaving the difference unchanged 
up to order a 3 . 

4 One might ask whether retaining the heavy quark mass in the kinematics (in the spirit of the "slow- 
rescaling" correction) might prove beneficial. The answer is that as long as the order of the parton distribution 
functions and the hard scattering coefficients are matched, this simply amounts to a shuffling of Mh/Q terms 
between the quark and gluon initiated terms, c/., Eq. ([|). 
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2 2.5 3 3.5 4 4.5 



Q (GeV) 

Figure 2: CTEQ4M charm quark density /h/ p (x, Q) and the approximate form fH/ P (x, Q), 
Eq. (|28|) , at x = 0.05 as function of Q- Since Eq. (|28|) is an order a] perturbative approx- 
imation while the evolution kernel for the CTEQ4M parton distribution includes order a 2 
terms, these curves do not match closely at threshold. 




2 2.5 3 3.5 4 4.5 



Q (GeV) 

Figure 3: CTEQ4M charm quark density /^/ p (x, Q) and an approximate form f^) p { x ^ Q) 
at x = 0.05 as function of Q. The approximate form is based on an analogue of Eq. (|28| ) 
in which appropriate order a 2 s terms are added. Since the calculation of fjjJ p (x,Q) uses 
the second order evolution kernel used for the CTEQ4M parton distributions, these curves 
match closely at threshold. 
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Q (GeV) Q (GeV) 

Figure 4: F^ for i = 0.1 as a function of Q as calculated using the ZM-VFN, FFN, 
ACOT, and S-ACOT prescriptions. The hard scattering coefficients are calculated to order 
ot\. Plot a) uses CTEQ4L parton densities and Plot b) uses CTEQ4M parton densities. 



quark has its proper mass but does not appear as a parton. For simplicity, we take jj, = Q. 
In our calculations, we evaluate the hard scattering coefficients F at order a s . Thus for the 
ZM-VFN, ACOT, and S-ACOT prescriptions, the terms displayed in Eq. fll5|) are included. 
The functions F are given by Eq. fl2|) for ZM-VFN, by Eqs. (0^171,0) for ACOT, and 
by Eq. (27) for S-ACOT. For the FFN prescription, there is only one term at order a s , as 
displayed in Eq. (ffijf). 

In Fig. 4 we show F£(x,Q) as a function of Q for x = 0.1. Then in Fig. 5 we show 
F£(x, Q) as a function of Q for x = 0.001. In each case we display results using both the 
CTEQ4L and CTEQ4M parton distributions. 

When we use the CTEQ4L parton distributions, we notice that there is a close match 
between the S-ACOT result and the FFN result near Q = m c . Based on the results of the 
previous subsection, we expect this matching to be degraded when we use CTEQ4M parton 
distributions because of the important role played by the order a 2 s term in the evolution 
kernel that is not matched in the lowest order calculation of the hard scattering function. 
This degradation is seen in the figures. 

In the asymptotic regime, Q 3> M, we find the S-ACOT result approximates the ZM- 
VFN result, as expected. 

We observe that the ACOT and S-ACOT prescriptions are effectively identical through- 
out the kinematic range. There is a slight difference in the threshold region, but this is small 
in comparison to the size of the /i- variation (not shown). Hence the difference between the 
ACOT and S-ACOT results is of no physical consequence. The fact that the ACOT and 
S-ACOT prescriptions match extremely well throughout the full kinematic range provides 
explicit numerical verification that the S-ACOT prescription fully contains the physics. 
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Figure 5: -F 2 C for x = 0.001 as a function of Q as calculated using the ZM-VFN, FFN, 
ACOT, and S-ACOT prescriptions. The hard scattering coefficients are calculated to order 
a\. Plot a) uses CTEQ4L parton densities and Plot b) uses CTEQ4M parton densities. 



7 Conclusions and outlook 

We have performed a numerical study of different prescriptions for dealing with the quark 
mass in heavy quark leptoproduction. We have seen that the simplest prescription f9j, S- 
ACOT, is numerically equivalent to the earlier ACOT prescription |4j]. 

The S-ACOT prescription is extensible order by order in a s . At 0{a\) we already 
find a significant simplification in the S-ACOT prescription as compared with the ACOT 
prescription. We expect that this simplification becomes even more dramatic at higher 
orders. 

We have not attempted to implement the S-ACOT prescription at order a 2 sl but we note 
that the NLO corrections in the the FFN prescription have been calculated 0, and that the 
leading (collinear) logarithms of the type a l s \og t (Q 2 /Mf I ) have been extracted in analytic 
form ||. Thus the S-ACOT subtraction term can be constructed as well. 

Finally, we emphasize that the choice of a prescription for dealing with quark masses 
in the hard scattering coefficients for deeply inelastic scattering is a separate issue from the 
choice of definition of the parton distribution functions. For all of the prescriptions discussed 
here, one uses the standard MS definition of parton distributions. 
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